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Abstract 

We carried out a first principles investigation on the microscopic properties of nickel-related 
defect centers in diamond. Several configurations, involving substitutional and interstitial nickel 
impurities, have been considered either in isolated configurations or forming complexes with other 
defects, such as vacancies and boron and nitrogen dopants. The results, in terms of spin, symmetry, 
and hyperfine fields, were compared with the available experimental data on electrically active 
centers in synthetic diamond. Several microscopic models, previously proposed to explain those 
data, have been confirmed by this investigation, while some models could be discarded. We also 
provided new insights on the microscopic structure of several of those centers. 

PACS numbers: 61.72.Bb, 71.55.-i,71.55.Cn 
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I. INTRODUCTION 



Diamond is a material which stands alone in nature, carrying a unique combination of elec- 
tronic, mechanical, thermal, and optical properties. Diamond is the hardest known natural 
material, having a large bulk modulus, high thermal conductivity and a large electronic band 
gap. Those properties make it a prototypical material to a number of applications, ranging 
from drilling and cutting tools to electronic devices to operate under extreme conditions 1 . 
More recently, new potential applications for doped diamond have been proposed, such as 
superconducting materials 2 and quantum computing^. 

There are two major methods with widespread use to grow macroscopic samples of syn- 
thetic diamond. Chemical vapor deposition (CVD) methods produce high quality diamond 
thin films grown over large areas. On the other hand, the high pressure-high temperature 
(HPHT) processes produce bulk diamond at relatively high growth rates and low costs. In 
those processes, samples are grown out of graphite, using 3d transition metal (TM) alloys 
(involving nickel, cobalt, and iron) as solvent-catalysts. Nickel is the only impurity that has 
been unambiguously identified in the resulting diamond. Such residual nickel impurities, 
either isolated or forming complexes with other defects, can generate several electrically 
and optically active centers^. Understanding the nature and microscopic structure of those 
centers is crucial in developing diamond-related technologies. Following the experimental 
identification of those centers, several microscopic models have been proposed to explain 
such data. However, there is still considerable controversy over a unified model which could 
explain most of those active centers in diamond. Here, we used first principles calculations 
to address this question. 

Over the last decade, nickel-related impurities in diamond have been investigated by 
several theoretical approaches^ - -. However, no investigation has provided a comprehensive 
picture of most of the nickel-related active centers identified in diamond so far. We used 
first principles total energy calculations, based on the full-potential linearized augmented 
plane wave methodology 10 , to investigate the structural and electronic properties of those 
centers in terms of the spin, symmetry, ground state multiplet, formation and transition 
energies, and hyperfine parameters. We focused our investigation on centers involving iso- 
lated nickel, either in interstitial or substitutional configurations, and complexes involving 
nickel and vacancies or dopants (boron and nitrogen). This paper is organized as follow: 
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in section [TT] we discuss the available electron paramagnetic resonance (EPR) experimental 
data of nickel-related defects in diamond. In section IHIl we present the methodology used 
in this investigation. Sections [TV] and [V] present and discuss the results in the context of 
experimental data. 

II. A SURVEY ON THE EXPERIMENTAL DATA 

Electron paramagnetic resonance and optical absorption measurements have been used 
to identify a number of nickel-related active centers in diamond, and have been recently 
reviewed^. The electronic properties of those centers have been analyzed in terms of either 
the Ludwig- Woodbury (LW) 11 or the vacancy^ 2 - model. According to the LW model, when 
a 3d n 4s 2 ion (1 < n < 9) occupies an interstitial site in a type- IV semiconductor, its 
4s electrons are transferred to the 3d orbitals, resulting in a 3d n+2 configuration. In the 
tetrahedral crystal field, the 3d states are split into e + t 2 irreducible representations. The 
threefold t 2 states lie lower in energy than the two-fold e states. This level ordering is the 
result of the octahedral crystal field, created by the next nearest neighbors of the impurity, 
which is stronger than the tetrahedral crystal field from the nearest neighbors. The same ion 
in a substitutional site would present a 3d 11-2 configuration, since four electrons are needed 
to bind with the four nearest neighboring host atoms. However, in this case the crystal 
field has tetrahedral symmetry, driving the e states to lie lower in energy than the t 2 ones. 
Additionally, in the LW model, the gap levels are filled according to the Hund's rule. The 
LW model is schematically presented in fig. [TJ 

For substitutional impurities, there is an alternative model, called the vacancy model, 
proposed by Watkins for TM elements near the end of the 3d, 4d, and 5d series^ 2 -. This 
model proposed that the electronic structure of the impurity resulted from a weak interaction 
between the impurity-related d-t 2 states and the t 2 vacancy-related ones, as represented 
in fig. [TJ The vacancy-related states came from the dangling bonds on the host atoms 
surrounding the vacant site into which the transition metal was inserted. As a result, the 
impurity band gap states would have a vacancy-like behavior. Although those two models 
were developed to describe the properties of 3d-transition metal impurities in silicon, they 
have been extensively used to explain the microscopic properties of those impurities in other 
semiconductors, such as nickel-related impurity centers in diamond. 
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Nickel in diamond has been detected in a tetrahedral symmetry with a spin S=3/2 by 
EPR 1 - and optical measurements^, and has been labeled W8 center. The microscopic 
model suggested for this center, based on either the LW or vacancy models, is an isolated 
substitutional nickel in the negative charge state (Ni~) in a 3d 7 configuration^. 

Two major active centers have been found in synthetic diamond, which have been associ- 
ated to interstitial nickel, labeled NIRIM-1 and NIRIM-2 centers^. The NIRIM-1 has been 
identified with a spin S=l/2 in a trigonal symmetry at low temperatures (T< 25 K), which 
switches to a tetrahedral symmetry at higher temperatures. This center was discussed in the 
context of the LW model, and interpreted as resulting from an isolated interstitial nickel in 
the positive charge state (Ni+ )S Since substitutional nickel in a positive charge state would 
give a spin S=5/2 according to the LW model, it was ruled out as a possible microscopic 
configuration for the NIRIM-1 center. More recently, independent investigations suggested 
that this center could, in fact, be formed by Ni+, giving a spin S=l/2&i£, indicating that 
the vacancy model is more suitable to describe this center. 

The NIRIM-2 center has been identified with a spin S=l/2 with a strong trigonal 
distortions^. The microscopic structure of this center is still the subject of controversy. It 
was initially associated to an interstitial nickel with an impurity or vacancy nearby^. More 
recently, this center has been proposed to be formed by a complex of nickel and boron^ or 
even by an isolated interstitial NA 

Post-growth annealing treatments introduce new active centers in the as-grown samples, 
which have been labeled NE centers^. It has been suggested that those NE centers involve 
nickel, nitrogen and vacancies. The NE4 center, that displays a symmetry and a spin 
S=l/2, has been tentatively associated to an interstitial Ni sitting in the middle position 
of a divacancy. This vacancy-nickel-vacancy unit (VNiV) would be aligned along a (111) 
direction. This configuration has also been labeled as NiCg in the literature, which rep- 
resents, besides the impurity, the six nearest neighboring carbon atoms. The LW modet^ 
was invoked to describe this configuration, in which the Ni impurity would donate six of its 
ten d-electrons to form bonds to the six neighboring carbon atoms. The remaining four 3d 
electrons should occupy a triplet d-orbital, according to Hund's rule. Since the center has a 
spin S=l/2, then this center should be in the negative charge state (V~NiV)~ , associated to 
a t\ electronic configuration. The NE4 center is the precursor to several other NE centers, 
which are formed by replacing nearest neighboring carbon atoms with nitrogen ones^, in 
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a NiC6- m N m (1 < m < 6) configuration. Recently, a center with rhombohedral symmetry 
and spin S=l, has been observed in diamond—. A (VNiV) configuration was suggested for 
this center, similar to that of NE4, but here in the neutral charge state. The NE1 center 
has a monoclinic symmetry and spin S=l/2, and has been suggested to be formed by a 
(VNiV) configuration plus two nearest neighboring nitrogen atoms^i. The NE8 center has 
a monoclinic symmetry and spin S=l/2, and has been suggested to be formed by a (VNiV) 
configuration plus four nearest neighboring nitrogen atoms 2 ^. 

EPR data has unambiguously shown that nickel can pair with boron and nitrogen impu- 
rities in diamond, forming new active centers. A center, labeled NOLI, has been identified 
with spin S=l and trigonal symmetry 2 ^. It has been suggested that this center is formed 
by an interstitial Ni 2+ (3d 8 ) impurity axially distorted by a boron (B~) along a (111) di- 
rection, with an unspecified interatomic distance between the impurities. A more recent 
examination of the trigonal boron-related NOLI center suggested a different model, that 
would involve substitutional nickel and boron, Ni+B°, with the acceptor boron in a next 
nearest neighboring site, with no covalent bonding between the impurities^. 

In samples with high concentrations of both nickel and nitrogen, other active centers 
have been identified. In addition to the NE centers, the AB5 center, with a spin S=l and 
trigonal symmetry, has been identified 2 ^. The microscopic model proposed for this center 
is a substitutional nickel (Nig - ) with a nearby substitutional nitrogen atom (N+). Table H] 
summarizes the properties of nickel-related EPR active centers in diamond, as well as the 
respective proposed microscopic models. 

Most of the microscopic models proposed in the previous paragraphs have been built 
based on an ionic model 24 , which has been proposed to describe the 3d-transition metal- 
acceptor pairs TMj -A~ in silicon^ (TM in a tetrahedral interstitial site in a positive charge 
state plus an acceptor A in a negative charge state). According to that model, the pair stable 
configuration corresponds to a classical system consisting of a TM+ electrostatically bound 
to a nearest neighbor A~ embedded in a dielectric medium 2 ^. Since the negatively charged 
acceptor has a closed shell, the electronic properties of the pair can be directly related to 
the positive TM ion placed in a screened Coulomb field. 
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III. METHODOLOGY 



We used the all-electron spin-polarized full-potential linearized augmented plane wave 
(FP-LAPW) method-, implemented in the WIEN2k package 26 . The calculations were per- 
formed within the framework of the density functional theory, using the Perdew-Burke- 
Ernzerhof exchange-correlation potential 27 . All the calculations were performed considering 
a 54-atom reference supercell. The methodology separates the crystalline space in two dis- 
tinct regions: the atomic and interstitial ones. The electronic wave functions were expanded 
in terms of spherical harmonics in the atomic regions and of plane waves in the interstitial 
ones. We chose all atomic spheres with a radius of R = 0.64 A. Therefore, 2R was much 
smaller than the crystalline interatomic distance of 1.54 A, such that atomic sphere overlap 
was avoided even in the case of large atomic relaxations. We used a 2 x 2 x 2 grid to sample 
the irreducible Brillouin zone as well as the T-point. 

Convergence in the total energy was tested by varying the number of plane waves de- 
scribing the electronic wave-functions in the interstitial region, a 7.0/R value provided con- 
verged results. Self-consistent interactions were performed until total energy and the to- 
tal charge in the atomic spheres changed by less than 10 -4 eV/atom and 10 -5 electronic 
charges/atom between two iterations, respectively. Additionally, the atomic positions were 
relaxed until the forces were smaller than 0.02 eV/A. All those approximations and conver- 
gence criteria have been shown to provide an accurate description of several defect centers 
in semiconductors 2 ^— . 

Formation and transition energies of all centers were computed using the procedure dis- 
cussed in Ref.— . The procedure required the total energies of the respective defect center 
and the chemical potentials of carbon, nitrogen and nickel. Those chemical potentials were 
computed using the total energy of carbon in a diamond lattice, nitrogen in a N 2 molecule, 
and nickel in a FCC lattice. In order to compute the hyperfme tensors, spin-orbit cou- 
pling was included in a second-variational procedure. Additional information concerning 
the calculation of hyperfme tensors is presented in Appendix. 
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IV. RESULTS 



We considered most of the proposed microscopic models for the electrically active centers 
described in table [U as well as, other possible models for those centers. Figure [2] represents 
the diamond lattice in the (110) plane, showing the possible sites in which the impurities 
could be placed in the beginning of each simulation. All atomic positions were later relaxed, 
according to convergence criteria discussed in the previous section. 

A. Isolated Nickel 

Substitutional nickel in diamond was considered in several charge states, with the results 
summarized in table [TTJ In the case of the neutral charge state (Nig), the center has no point- 
symmetry (Ci), and presents a spin S=l. This configuration is only 0.1 eV more stable than 
the center in a symmetry (Nig*). Figure [3] displays the induced energy eigenvalues in the 
gap region for substitutional nickel impurity. The gap states of Ni s are vacancy-like orbitals, 
consistent with the vacancy model^. 

Table [III also presents the results for interstitial nickel. In the positively charge state 
(Ni^), the center was initially simulated in a trigonal (C3„) symmetry, in order to check if 
that configuration could explain the properties of the NIRIM-2 center™. In that symmetry, 
it presented an effective spin S=l/2 and an 2 E multiplet ground state. By releasing the 
symmetry constraint, there was an energy gain of about 0.2 eV, and the center distorted 
to a Cih symmetry. This symmetry lowering was very small, corresponding to a distortion 
on the nickel atom of only 0.06 A toward one of its second nearest neighbors, breaking the 
trigonal symmetry. Figure H] compares the electronic structure of the Ni^ center in both 
symmetries, showing that although the symmetry lowering was small, there were strong 
effects in the electronic structure of the center. These results show that the electronic 
structure of interstitial nickel cannot be described by the LW model^, since the 3d-nickel 
related states remain resonant in the valence band, leaving a hole in the perturbed valence 
band top. In trigonal symmetry, the valence band top of the diamond crystal splits into an 
ai state, resonant in the valence band, and an e state, occupied by three electrons, inside 
the gap. In the C\h symmetry, the e gap states split further in an a' and an a" ones. 
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B. Ni- vacancy complexes 

We initially considered an interstitial nickel paired with a nearest neighboring vacancy 
(NiiV), as suggested as a stable configuration in several experiments^. However, this config- 
uration was unstable and the impurity moved toward the vacant site, forming a substitutional 
nicked. We additionally considered a substitutional nickel paired with a nearest neighboring 
vacancy (Ni s V), but in the final relaxed structure, the nickel remained in the middle position 
between two vacancies (VNiV). Figure |5] presents the induced energy eigenvalues of this 
last complex and table [II] presents the respective properties. 

The electronic structure of the V^NiV^ complex cannot be described by the LW model, as 
it has been recently suggested^. Our results indicate that the relevant electronic properties 
of this center should be associated to divacancy-like orbitals, which appeared in the gap, 
while the Ni-related orbitals remained resonant and inert inside the valence band. On the 
other hand, the electronic structure is well described by the crystal field theory, in which the 
electronic states can be interpreted as resulting from an interaction between the divacancy 
states and those of the Ni atom. The one-electron ground state structure of a diamond 
divacancy in D 3 ^ symmetry has the «2u a ig e u e g configuration. In that symmetry, the Ni 3d 
energy levels are split into 2e g + ai g . When a Ni atom is placed in the middle position of 
a divacancy, its e g energy level interacts with the carbon dangling bonds, leaving a fully 
occupied non-bonding t 2g -like (e g + ai g ) orbital inside the valence band. On the other hand, 
the Ni e g state interacts with the divacancy e g gap level, leaving the e g -bonding level in the 
valence band and the e g -anti-bonding one unoccupied in the gap. The relevant electronic 
properties of this center are related to the e u divacancy-like orbital, which remained in the 
gap bottom. In the positive and negative charge states, the symmetry lowering (D^d — > C2h) 
is very weak and the splitting in the e u -related states is smaller than 0.1 eV. 

C. Ni-B complexes 

We now consider complexes involving nickel and substitutional boron, which could poten- 
tially lead to a trigonal symmetry, to be consistent with proposed models for Ni-B centers 
presented in table [I] For interstitial nickel-substitutional boron pairs, we considered three 
microscopic configurations, according to figure EJ NijB s pair, with Ni and B respectively 
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in sites 6 and 1; Ni;CB s pair, with Ni and B respectively in sites 6 and 4; and Ni;(g)B s 
with Ni and B respectively in sites 7 and 1. Table II I II presents the properties of those 
centers in several charge states. The (Ni;B s ) * complex, has a degenerate configuration in 
C^ v symmetry, coming from the partially occupied e state, favoring a symmetry lowering 
to Ci. The distance between the Ni and B in the pairs is crucial for the final properties of 
those centers, as evidenced by the electronic structure of those three centers, in the same 
charge state, shown in figure [6) The major difference emerges on the character of the highest 
occupied level in the center. While for the NijB s and NiiCB s pairs, this level has a localized 
Ni 3d- related character, for the Nii®B s , this level is essentially delocalized. For this last 
center, the distance between the impurities is so large that the center can be well described 
by an ionic model, in which the role of boron is only to accept an electron from the nickel 
impurity. Therefore, the electronic structure of this complex can be well described as an 
isolated interstitial nickel in 2+ charge state (Ni 2+ ). 

For substitutional nickel-substitutional boron, we considered two structural configura- 
tions, according to figure EJ Ni s B s pair, with Ni and B respectively in sites 3 and 4; Ni s ®B s 
pair, with Ni and B respectively in sites 1 and 5. Figure [7] presents the energy eigenvalues 
of those pairs and table [TTT1 summarizes their properties. In the Ni s ®B s centers, boron is far 
from the nickel impurity, working as just an acceptor, such that the electronic configuration 
resembles that of isolated substitutional Ni impurity, shown in figure [3j For the Ni s B s pair, 
boron plays a more important role, affecting the electronic structure of the center, although 
the magnetic properties of this center are associated with partially occupied energy levels 
with prevailing nickel character. 

D. Ni-N complexes 

The nickel and nitrogen complexes in diamond are generally formed as result of high- 
temperature thermal annealing, in which nitrogen impurities become highly mobile and end 
up pairing with the less mobile nickel ones. We considered centers with nickel in interstitial, 
substitutional, and divacancy sites complexing with nitrogen. According to fig. [21 the NiiN s 
center has Ni and N atoms respectively in sites 6 and 1 and the Ni s N s center has Ni and N 
atoms respectively in sites 3 and 4. Table [IV] presents the results for Ni-N pairs. 

Figure M describes the electronic structure of Ni-N complexes as resulting from an inter- 
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action of the energy eigenvalues of the impurities in isolated configurations. For the (Ni s N s )° 
complex, there is a hybridization between 3d-related gap states of Ni s with those 2p-related 
ones of N s . This indicates a typical covalent interaction between those two impurities^i^. 
Despite this hybridization, the highest occupied energy level in the (Ni s N s )° complex has 
an e representation, with a prevailing 3d character. On the other hand, for the (Ni;N s ) 
complex, the electronic structure results from a weaker interaction between the states of 
the isolated impurities, more consistent with an ionic model^. In this last case, the highest 
occupied level has both 3d-related Ni and 2p-related N characters. 

Table IIVI presents the results for complexes involving nickel and nitrogen impurities in 
a divacancy site. We considered two possible configurations, according to the proposed 
models for the NE1 and NE8 active centers described in table HI The N VNi VN complex 
involves the precursor FNiF plus two nitrogen atoms in diametrically opposed positions, 
replacing two of the nickel six nearest neighboring carbon atoms. The N2 VNi VN2 complex 
has four substitutional nitrogen atoms, replacing four of those nearest neighboring carbon 
atoms. The electronic structure of those two centers shows a strong covalent interaction 
between the divacancy-related orbitals and the nitrogen-related ones, which is similar to 
what is observed for complexes involving cobalt-nitrogen complexes in diamond—. Nitrogen 
incorporation into the precursor substantially alters the electronic structure of that center. 
This shows that the current interpretation, in which nitrogen atoms play a role of only 
donating electrons to the precursor, is not valid. 



V. DISCUSSIONS AND CONCLUSIONS 



The results are now discussed in the context of the experimentally identified active cen- 
ters observed in synthetic diamond. In a previous investigation, we have shown that the 
microscopic configuration of a substitutional nickel in the negative charge state (Ni~) was 
consistent with the properties of the W8 center-^, including symmetry, spin, and hyperfine 
parameters. Additionally, we have shown that the previously proposed microscopic models 
for the NIRIM centers, described in table HI based on interstitial nickel impurities, were 
not consistent. For example, the NIRIM- 1 center could be better explained by an isolated 
substitutional nickel in the positive charge state (Ni+)&^, although this investigation finds 
a number of other configurations which are also consistent with the NIRIM-1 symmetry and 
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spin. 

For the NIRIM-2 center, a direct comparison between theory and experiment was more 
complicated. One of the proposed microscopic models for the NIRIM-2 center was an inter- 
stitial nickel with a nearby vacancy^. Theoretical investigations showed that this configura- 
tion is unstable 8 -, such that the interstitial nickel migrates toward the vacant site, becoming 
a substitutional impurity. This would be fully expected considering defect energetics, since 
the formation energy of substitutional nickel is considerably lower than that of an intersti- 
tial one. We have previously suggested that isolated interstitial nickel in the positive charge 
state (Ni+) could explain some of the properties of the NIRIM-2 center 8 . This investigation 
shows that Ni^ is unstable in trigonal symmetry, lowering to a Ci one. However, the energy 
gain from this symmetry lowering is only 0.2 eV, and the final configuration is not far from a 
trigonal symmetry. Recently, it has been proposed that NIRIM-2 should involve interstitial 
nickel with a next nearest neighboring boron atom^. The results for this proposed configu- 
ration, (NiiCB s )°, are fully consistent with the experimental data for NIRIM-2 in terms of 
symmetry and spin. Another center, involving boron and substitutional nickel, (Ni s ®B s )°, 
also provides results consistent with experimental data of NIRIM-2. This last configuration 
would be a strong candidate to explain the NIRIM-2 center since it involves substitutional 
nickel, and formation energy is considerably smaller than that for a pair involving intersti- 
tial nickel. A definite answer on the NIRIM-2 microscopic model could be achieved if future 
experiments could resolve the Ni-hyperfine parameters, since according to table IHIl those 
parameters are considerably different for those two configurations. 

Another center has been associated to nickel-boron pairs. The NOLI center, probably 
the same as the NIRIM-5 center, has been found in heavily boron-doped diamond^ 1 ^. The 
center has trigonal symmetry and S=l. By inspection of our results, the (Ni;B s ) + complex, 
suggested by^ as the microscopic structure of this active center, is fully consistent with the 
experimental data. Another complex involving interstitial nickel, (NigB s ) + , is also consistent 
with experimental data. Complexes involving substitutional nickel could also describe the 
properties of the NOLI center. The (Ni s B s ) + complex, suggested in Ref.— , is diamagnetic 
and cannot explain the NOLI results. However, the same complex in a negative charge 
state, (Ni s B s )~, is fully consistent with the experimental data. Although this center has been 
only observed in heavily boron-doped diamond, our results indicate that nickel in isolated 
configurations, Ni^ - or Nif + , are also consistent with the experimental data. In the case of 
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the Nij , the high concentration of boron would only be required to place the Fermi level 
near the valence band top to access the 2+ charge state, and not necessarily participating 
into the complex. In order to clarify this, EPR experiments should be performed to observe 
the hyperfine parameters in nickel and boron nuclei. 

For the NE4 centers, experiments^ have suggested a microscopic structure given by a 
nickel impurity in a divacancy site in the negative charge state (VNiV)~ . Our results for this 
configuration give a trigonal symmetry and a spin S=l/2, both results consistent with the 
experimental findings. The NE4* center—, from table HJ has been suggested to be formed 
by (VNiK) . Our results corroborate that suggestion, although they indicate a trigonal 
symmetry, while experiments suggested a rhombohedral one. Moreover, our calculations 
found a hyperfine parameter (Aj_) of 60 MHz in the nearest neighboring carbon atoms, very 
close to the experimental value of 79 MHz 20 . 

The NE1 and NE8 centers have been suggested to be formed by nickel-nitrogen complexes 
in a divacancy site 2 ^. Our results, in terms of spin and symmetry, for the (NFNiFN) - 
and (N 2 FNiFN 2 ) + complexes are fully consistent with the experimental data and the 
proposed microscopic configurations. However, according to table WV\ the (NFNiFN) 4 " 
and (N2FNiV r N2)~ complexes also provide results consistent with those data. However, 
experiments 21 could not resolve the Ni-related hyperfine fields, in order to compare with the 
values presented in table IIVI On the other hand, those experiments have identified hyper- 
fine fields in the nitrogen and the nearest neighboring carbon nuclei. For the NE1 center, 
the experimental values for those fields are Ay (N) = 59, Aj_(N) = 40, A||(C) = 49, and 
A ± (C) = 31 MHz. For the (N VNi VN)~ complex, our results provide Ay(N) = 42 and 
Aj_(N) = 17 MHz and negligible values in the carbon nuclei. For the (NFNil/N) + complex, 
our results provide A||(C) = 92 and Aj_(C) = 40 MHz and negligible values in the nitrogen 
nuclei. Therefore, it was not possible to make a final remark on the microscopic structure 
of the NE1 center. For the (N2VN1VN2) 4 " complex, hyperfine fields in the nitrogen and 
carbon nuclei are fully consistent with experimental values of the NE8 center. Finally, the 
(Ni s N s )~ complex has been proposed as the microscopic structure of the AB5 center—. From 
all the complexes involving nickel and nitrogen considered here, that configuration was the 
only one consistent with the experimental results of the AB5 center. 

In summary, we have performed a theoretical investigation on nickel-related complexes 
in diamond, in terms of electronic structure and hyperfine fields. We have explored several 
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microscopic configurations that could explain the experimental data on EPR active centers 
in synthetic diamond, confirming or discarding some of the previously proposed microscopic 
models and suggesting new ones. These results provide a comprehensive picture on Ni- 
related active centers in diamond using a single theoretical methodology 

Acknowledgments 

The authors acknowledge partial support from Brazilian agency CNPq. The calculations 
were performed using the computational facilities of the CENAPAD and the LCCA-CCE 
(Universidade de Sao Paulo). 

APPENDIX: CALCULATION OF THE HYPERFINE TENSORS 

The EPR data can provide important information related to electrically active centers 
in semiconductors, such as symmetry, spin, gyromagnetic factor, and in some stances the 
atomic composition of those centers. The hyperfine spectrum of a center results from an 
interaction between nuclear magnetic moments (/ij) and the moments of unpaired electrons. 
The crystalline field, in which the impurity (or other defects) is immersed, is generally enough 
strong to quench the respective orbital moment. However, it has been shown that in the case 
of transition metal impurities, the orbital moment is not fully quenched by the crystal field, 
generating, in some cases, a large energy anisotropy. Most of the theoretical investigations 
have neglected this anisotropic contribution, but it is very important for systems such as 
those investigated here. 

The hyperfine fields were computed using the implementation from the WIEN2k 
package^ that uses a scalar-relativistic approximation™ According to that approach, the 
hyperfine magnetic field (EW) is computed considering three components: the Fermi contact 
(B c ), the dipolar (Bd ip ), and the orbital (B orb ) terms. 

Bhf = B c + B dip + B orb , (1) 

These three components are given in terms of the the angular (L) and spin (S) electronic 
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moments (in h unities) and the Bohr magneton (j3 e = eh /2m): 



B dip = -^c/3 c ($ 



sir) 



S-3(S-r] - 



B 



orb 



si?) 



L 



(2) 
(3) 

(4) 



where <£> is the relativistic large component of the wavefunction and S(r) is the reciprocal 
relativistic mass enhancement: 

i + e ~ v ^y 1 



Sir) 



2mc 2 



(5) 



where e and V(r) are respectively the kinetic energy and the Coulomb potential. 
rh av is the average nuclear magnetization, 

= J 5 T (r')m(r')dr' = 

5 T (r')($ \ aS{f-r') | $ ) dr ' 



where 5^{r') is given in terms of the Thomas radius (pp = Ze 2 /mc 2 ): 

Sri?) 1 



(6) 



(7) 



Arrr 2 [2r(l + e/2mc 2 ) + r T ] 
and c? are the Pauli matrices. 

The splitting in energy resulting from the interaction between the hyperfme magnetic 

— * 

fields (Bhf) and fl\ is described by: 



E = — p,j ■ B hf . 



(8) 



This splitting in energy may be described in terms of a spin Hamiltonian (H). The eigen- 
values of this Hamiltonian provide information on the separation between absorption lines 
in the magnetic spectra. 

H = J • A • I = (L + S) • A • I (9) 

where I is the nuclear spin and A is a (3 x 3) hyperfme interaction tensor. 
The hyperfme interaction tensor has the following components Ay: 

a^y + 4 ip + a°/ b , with (10) 



A 



£a? p = and 
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In an experiment, when the direction of the external static magnetic field (h = 
sin 9 cos ipi + sin^sin^j + cos 9 k) is varied with respect to the sample axis, the relevant 
quantity is the projection of the hyperfine interaction tensor in that direction: 

A(9, <p) = n • A • h = 
= An sin 2 9 cos 2 ip + (A i2 + A 2 i) sin 6 cos ip sin ip + 
+ A 22 sin 2 9 sin 2 ip + (A 2 3 + A 32 ) cos 9 sin 9 sin ip + 
+ A 33 cos 2 + (Ai 3 + A 3i ) cos # sin # cosy?. (11) 

By choosing a convenient set of six directions, i.e. six sets of (9, ip), the values of A(9, <p) in 
those directions allow to build the hyperfine interaction tensor. It can be later diagonalized 
to obtain the three principal values, also called hyperfine parameters (Ai, A 2 , and A 3 ), and 
their respective eigenvectors. 

•O- 

The hyperfine tensor A is given in terms of: 

A = a c l+B + C, (12) 

where 1 is the unitary tensor, and a c is the contact term, B is a traceless anisotropic tensor 
related to the dipolar interaction, and C is an anisotropic tensor related to the orbital 
interaction. 

If the angular magnetic moment is quenched, the isotropic part of the hyperfine tensor 
is exactly the Fermi contact interaction and the anisotropic part is the dipolar interaction. 
However, if the angular moment is not quenched, there should be a contribution from this 
interaction to the hyperfine tensor changing both the dipolar and contact terms. In this 
investigation, we observed that the hyperfine orbital field is generally relevant and cannot 
be neglected. This was result of the spin-orbit coupling in the 3d localized orbitals, which 
are deformed due to the crystalline field. 
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TABLE I: Experimental data on the electrically active centers of Ni-related defects in diamond. 
The table presents the symmetry, spin, and proposed microscopic model. X represents an unknown 
specie (a vacancy or impurity), and V represents a vacancy. A list of additional relevant active 
centers in diamond can be found in Reference^. 



Label 


Sym. 


S 


Model 


W8 


tetrahedral 


3/2 


Ni- (°) 


NIRIM-1 


trigonal (T< 25K) 


1/2 


Ni+ (*), Ni+ W'M 


NIRIM-2 


trigonal 


1/2 


Ni+-X W, Ni+CB S - M) 


NE4 


trigonal 


1/2 


(VNiV) - (e) 


NE4* 


rhombohedral 


1 


(vmv)° (/) 


NE1 


monoclinic 


1/2 


(WNiVN) - fo) 


NE8 


monoclinic 


1/2 


(N 2 yWN 2 ) + (9) 


NOLI 


trigonal 


1 


Ni+Bg ( c ), Nif + B s - (a) 


AB5 


trigonal 


1 


Ni h 2 -N s + 00 



( a ) Reference^, ^ Reference^, ( c ) Reference^, ^ Reference^, 
^ Reference^, Reference^, ^ Reference^, CO Reference^. 
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TABLE II: Results for isolated Ni and Ni-divacancy complexes in diamond: symmetry, spin (S), 
multiplet ground state ( 2S+1 T), formation energies (Ep), and transition energies (E t with relation 
to the valence band top e v ). Here ep is the Fermi energy. The table also presents the calculated 
hyperfine parameters (Ai, i = 1, 2, 3) in the 61 Ni nucleus. Energies and hyperfine parameters are 
given in eV and MHz, respectively. 



Center 


Sym. 


S 


2S+lp 


Ep 


E t 


Ai 


A 2 


A 3 


M s 2+ 


T d 







3.9 + 2e F 


2.0 (2 + /+) 


— 


— 


— 


Ni+ 


Cs v 


1/2 


2 Ai 


5.9 + €F 


2.6 (+/0) 


123 


-36 


-36 


Ni°* 

S 


^ ou 


1 


3 E 


8.6 




52 


9 


9 


Ni° 

S 


Ci 


1 


3 A 


8.5 




45 


18 


4 


Nir 

s 


T H 

a 


3/2 


4 A 2 


11.5 - e F 


3.0 (0/-) 


18 


18 


18 


Ni 2 " 




1 


3 A 


15 5 - 2ep 


4 o (-12-) 


-99 


21 


21 


Ni 2 + 




1 


3 A 


15.5 + 2e F 


0.6(2 + /+) 


32 


2 


2 


Ni+* 




1/2 


2 E 


16.3 + e F 




29 


15 


15 


Ni+ 


Cih 


1/2 


2 A 


16.1 + e F 


1.1 (+/0) 


66 


19 


17 


Ni? 


T d 





1 A 1 


17.2 










(vmv) + 


C 2 /i 


1/2 


2 A 


5.2 + e F 


0.2 (+/0) 


51 


17 


14 


(vmv)° 




1 


3 A 2u 


5.4 




6 


33 


33 


(vmv)~ 


C 2 h 


1/2 


2 A 


6.2 - e F 


0.8 (0/-) 


18 


-52 


-22 


{vmv) 2 - 


D 3d 





X A 


7.3 - 2e F 


1.1 (-/2-) 
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TABLE III: Results for nickel-boron complexes in diamond: symmetry, spin (S), multiplet ground 
state ( 2S+1 r), formation (E F ) and transition energies (E t with relation to e v ). The table also 
presents the calculated hyperfme parameters (A;) in the 61 Ni nucleus. 





Sym. 





2S+l-n 

1 


H;p 








An 


(Ni;B s ) + 




1 


3 A 


14.8 + e F 


1.1 (+/0) 


7 


35 


35 


(NiiB s )°* 


C% v 


1/2 


2 E 


16.1 




76 


-62 


-62 


(NiiB s )° 


Ci 


1/2 


2 A 


15.9 




40 


-56 


-27 


(NiiB s )- 


C% v 





1 A 1 


17.3 - e F 


1.4 (0/-) 


- 


- 


- 


(NiiCB s )+ 







1 A 1 


15.4 + e F 


1.8 (+/0) 


— 


— 




(NiiCB s )° 




1/2 


2 Ai 


17.2 




21 


-14 


1 ; 


(NiiCB s )- 







1 A 1 


18.2 - e F 


1-4 (+/-) 


- 


- 


- 


(Nii0B s )+ 




1 


3 A 


15.9 +e F 


0.5 (+/0) 


7 


-17 


-17 


(Nii®B s )° 


Ci 


1/2 


2 A 


16.4 




-5 


50 


46 


(Nii®B s )- 


C3V 





1 A 1 


17.9 - e F 


1.5 (0/-) 


- 


- 


- 


(Ni s B s )+ 


C3V 





1 A 1 


3.5 + e F 


2.6 (+/0) 








(Ni s B s )° 


Ci 


1/2 


2 A 


6.1 




-99 


55 


-32 


(Ni s B s )- 


C3V 


1 


3 A: 


9.0 - e F 


2.9 (0/-) 


85 


1 


1 


(Ni s ®B s )+ 


C3V 





1 A 1 


4.7 + e F 


2.1 (+/0) 








(Ni s ^B s )° 


C3 V 


1/2 


2 A 1 


6.8 




-120 


48 


48 


(Ni s ®B s )" 


Ci 


1 


3 A 


9.2 - e F 


2.4 (0/-) 


42 


18 


8 
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TABLE IV: Results for nickel-nitrogen complexes in diamond: symmetry, spin (S), multiplet ground 
state ( 2S+1 r), formation (Ep) and transition energies (E t with relation to e v ). The table also 
presents the calculated hyperfine parameters (A;) in the 61 Ni nucleus. 



Center 


Sym. 


s 


9S-I-1 t~i 


Ep 


E t 


Ai 


A 2 


A 3 


CNi-N ) + 


^Sv 


o 


X At 
■"■1 












CNi-N "1° 




1/2 


2 At 
■"■1 


19.0 




82 


14 


14 


(NiiN s )- 







X A X 


22.5 - e F 


3.5 (0/-) 








fNicN^" 1 " 




o 


x Ai 


6.0 + £f 


3 1 f + /q) 








CNi N "1° 


Hi 


1 /2 


2 A 


q 1 

C/. -L 




-i i n 


7D 


-3D 

OU 


CNi N "l - 




i 

j_ 


3 Ai 




,5.0 I^U/ J 


oo 


9 





(NFNiFN)+ 


C 2 /i 


1/2 


2 A 


3.2 + e F 


1.3 (+/0) 


28 


-20 


-24 


(NFNiFN) 


C 2 fe 





X A 


4.5 










(NFNiFN)- 


C 2 fe 


1/2 


2 A 


7.5 - e F 


3.0 (0/-) 


177 


58 


43 


(N 2 FNiFN 2 )+ 


C 2 ft 


1/2 


2 A 


2.7 + e F 


3.7 (+/0) 


-167 


-18 


-15 


(N 2 FNiI/N 2 )° 


C 2 fe 





X A 


6.4 










(N 2 FNiFN 2 )- 


C 2 fe 


1/2 


2 A 


10.7 - e F 


4.3 (0/-) 


1 


-5 


-4 
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LW 


Watkins 


Nii Ni s 


Ni 


Ni s V 


c=fc= u =f= 


3d(e) / 






/ ' 
i ' 

. i 
i 


3d 10 3d 6 


i 


3d(t 2 ) v 


i 



FIG. 1: Schematic representation of the gap states for an isolated interstitial (Nii) or substitutional 
(Ni s ) nickel in neutral charge state, according to LW— and Watkins (vacancy)^ models. The f 
and 1 arrows represent the spin up and down, respectively. Gray regions represent the valence 
and conduction host bands. For simplicity, the system is considered in a tetrahedral symmetry, 
neglecting distortions. 
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FIG. 2: Representation of the diamond lattice in the (110) plane. Carbon atoms are represented 
by black circles. Labels, from 1 to 5, indicate the crystal site positions where the impurities could 
be placed. The figure also shows two tetrahedral interstitial sites (6 and 7), represented by the 
symbol, in the [111] direction. 
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Ni s + Ni s + Ni s Ni s Ni s 2 



t e t 

l 2 (i5) a lt 



(a) (b) (c) (d) (e) 

FIG. 3: The Kohn-Sham spin-polarized energy eigenvalues (around the T point) representing the 
3(i-related Ni levels in the gap region for isolated substitutional nickel in different charge states: 
(a) Nig + , (b) Ni+, (c) Nig, (d) Ni~, and (e) Nig~. Levels with spin up and down are represented 
by t and 4- arrows, respectively. The occupation of the gap levels is given by the number of filled 
circles. Numbers in parenthesis represent the ci-character percentage of charge inside the Ni atomic 
sphere. 
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FIG. 4: The energy eigenvalues representing the 3<i-related Ni levels for isolated interstitial nickel in 
the positive charge state in Cs v and Ci^ symmetries. Levels with spin up and down are represented 
by t and 1 arrows, respectively. The occupation of the gap levels is given by the number of filled 
circles. The numbers in parenthesis represent the (i-character percentage of charge inside the Ni 
atomic sphere. 
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<yNiV) + 


(VNiV) 


(VNiV)~ 


(VNiV) 2 ' 


leV 






(21) 


eg* W 

e gt (24) 


(17) ... 


(18) 


■y' 


" (25) 


" (18) 






FIG. 5: The energy eigenvalues representing the 3d-related Ni levels in the gap region for the Ni- 
divacancy complexes: (a) (VTW)+ (b) (VNiV) , (c) (VNiV)-, and (d) (VNiV) 2 ~ centers. The 
occupation of the gap levels is given by the number of filled circles. The numbers in parenthesis 
represent the d-character percentage of charge inside the nickel atomic sphere. Levels with spin up 
and down are represented by f arid 1 arrows, respectively. 
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(Niftf (Nii®B s ) + (NiiCB s ) + 

leV 

(7) 

a (22) 

e l \ (6) 

e mm ^ •• *> m i 

(12) * mm - (7) 

(30) 

(a) (b) (c) 

FIG. 6: The energy eigenvalues in the gap region for three configurations involving interstitial 
nickel-substitutional boron complexes in the positive charge state: (a) (Ni;B s ) + , (b) (Nii(g)B s ) + , 
and (c) (NiiCB s ) + . The occupation of the levels is given by the number of filled circles. The 
numbers in parenthesis represent the d-character percentage of charge inside the nickel atomic 
sphere. 
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FIG. 7: The energy eigenvalues for the substitutional nickel-substitutional boron complexes in two 
configurations: (a,b,c) for the Ni s B s complex and (d,e,f) for the Ni s (g) B s complex. The occupation 
of the gap levels is given by the number of fillecl^circles. The numbers in parenthesis represent the 
d-character (o-character) nercentaee of charge inside the nickel (boron) atomic sohere. Levels with 
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FIG. 8: The energy eigenvalues for the (Ni; N s ) and (Ni s N s ) complexes in the neutral charge state. 
The figure shows that the electronic structure of those centers results from hybridization between 
the 2p nitrogen with 3d nickel levels coming from its precursors in isolated configurations. 
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